Current-Induced Dynamics in Almost Symmetric Magnetic Nanopillars 
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Magnetic nanodevices usually include a free layer whose configuration can be changed by spin- 
polarized current via the spin transfer effect, and a fixed reference layer. Here, we demonstrate 
that the roles of the free and the reference layers interchange over a small range of their relative 
thicknesses. Precession of both layers can be induced by spin transfer in symmetric devices, but the 
dynamics of one of the layers is rapidly suppressed in asymmetric devices. We interpret our results 
in terms of the dynamical coupling between magnetic layers due to spin transfer. 

PACS numbers: 85.75.-d, 75.60.Jk, 75.70.Cn 



Spin torque (ST) exerted by spin-polarized current 
I on a nanopatterned magnetic layer Fi can change 
its magnetic configuration or induce magnetic dynam- 
ics [l|, 0]. The polarization of I is induced by another 
magnetic layer F2, which is usually made much larger 
than Fi to minimize the effects of ST on this layer. How- 
ever, it is presently not known what makes a specific 
magnetic layer behave as a fixed polarizer or a free layer 
driven by ST. Significant effects of ST on both layers 
can be expected if their dimensions are similar, but the 
current-induced behaviors of such symmetric structures 
are not well understood. Magnetoresistance (MR) mea- 
surements of symmetric nanopillars indicated dynamics 
for both polarities of /, interpreted as independent exci- 
tations of each layer by the appropriate direction of I [3j • 
However, calculations Q and measurements of current- 
induced reversal statistics in symmetric nanopillars |5j 
suggest that ST-induced coupling results in simultane- 
ous incoherent dynamics of both layers. 

Significant effects of ST on the polarizing layer are also 
possible in asymmetric nanopillars [6(. Narrow spectro- 
scopic peaks at large / were attributed to the dynam- 
ics of the polarizer [7|, but their precise nature and the 
mechanism of excitation are not known. Spectroscopy 
of current-induced dynamics in nanopillars with similar 
dimensions of Fi and F2 can potentially clarify these ef- 
fects. Here, we report spectroscopic measurements of 
current-induced dynamics in almost symmetric nanopil- 
lars, where decoherence caused by ST-induced coupling 
between magnetic layers was reduced due to their dif- 
ferent dynamical properties. We demonstrate coher- 
ent precession for both polarities of / in symmetric de- 
vices. However, even slightly asymmetric devices exhibit 
a rapid suppression of precession for one of the current 
polarities, indicating significant effects of dynamical cou- 
pling between layers. 

Multilayers Cu(60)FiCu(8)F 2 Au(10)Cu(100), where 
thicknesses arc in nm, were fabricated on oxidized silicon 
by a technique described elsewhere 0. Only F 2 =Py(5) 
(Py=Ni 80 Fe2o) and part of Cu(8) spacer were patterned 
into an elliptical shape with approximate dimensions of 
100 nm x 50 nm, while ¥ i=Py {tp y ) was left extended 



with dimensions of several /j,m. This sample geometry 
was verified by MR measurements, which showed neg- 
ligible coupling between Fi and F 2 . All measurements 
were performed at room temperature (RT). The sam- 
ples were contacted by coaxial microwave probes, which 
were connected through a bias tee to a current source, 
a lock-in amplifier, and a spectrum analyzer through a 
broadband amplifier. The microwave signals were ad- 
justed for the frequency-dependent gain of the amplifier 
and losses in the cables and probes, determined with a 
microwave generator and a power meter. / > flowed 
from Fi to F 2 . The in- plane field H was rotated by 40° 
with respect to the nanopillar easy axis, unless specified 
otherwise. The rotation enabled spectroscopic detection 
of precession, and had an important effect on precession 
coherence Q. 

The effects of ST on Fi relative to F 2 were varied by 
changing its thickness tp y . We discuss current-induced 
behaviors of three samples with tp y = 2,3, and 4 nm, 
labeled Py2, Py3, and PyA. We shall see below that, 
because of the extended geometry of Fi, Py3 behaves as 
a symmetric sample with respect to the effects of cur- 
rent on Fi and F 2 , while the behaviors of Py2 and PyA 
exhibit opposite asymmetries. The difference Ai? be- 
tween resistance Rap in the antiparallel (AP) configura- 
tion and resistance Rp in the parallel (P) configuration 
was 0.082 fi, 0.124 O, and 0.134 Q, for Py2, Py3, and PyA 
respectively, consistent with stronger polarizing proper- 
ties of increasingly thicker Fi. The behaviors reported 
below were verified for at least two devices of each type. 

The power spectral density (PSD) data for sample Py2 
exhibit two harmonically related spectral peaks caused 
by current-induced precession at / < (Fig. QJa)). The 
onset of the microwave peaks above the measurement 
noise floor of 1 fW /MHz up to 8 GHz coincides with 
a small peak in differential resistance dV/dl at Iq — 
-2.3 mA, Fig. [T^b). A step in dV/dl at J+ = 3.8 mA 
corresponds to the onset of 1// noise at I > 0, char- 
acteristic of thermal transitions between different stable 
magnetic configurations of the bilayer. 

In contrast to Py2, the spectroscopic data for sam- 
ple Py3 exhibit precession peaks for both polarities of I 
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FIG. 1: (Color online) (a) PSD vs frequency / and J, and (b) 
5R = dV/dI(I) - dV/dI(0) for sample Py2 at H = 600 Oe. 
(c),(d) same as (a),(b), for Py3. (e),(f) same as (a),(b), for 
PyA, (g),(h) same as (e),(f), at H = 400 Oe. 



(Fig. [He)). Their onsets coincide with bumps in dV/dl 
at I c = -A.A mA and l£ = 3.8 mA (Fig. fQd))- The 
smallest full width at half maximum (FWHM) of the 
first harmonic at I < is 15 MHz, similar to that seen 
in point contacts [10|, but significantly smaller than in 
nanopatterned layers at RT [8(. The first harmonic at 
I > has the smallest FWHM of 164 MHz, and is super- 
imposed on 1// noise similar to that seen in Py2. The 
frequency of the first harmonic at the onset current Iq 
is lower than the precession frequency at Iq. Based on 
these differences between the dynamical features at / < 
and I > 0, we tentatively attribute the peaks at / < 
to precession of the extended Fi , and the peaks at I > 
to precession of the nanopatterned F2 . We generally ex- 
pect that the current-induced dynamics of F2 is localized 
in the area somewhat exceeding the dimensions of the 
point contact. The larger precession frequency of Fi can 
be explained by the exchange interaction between the 
precessing part of Fi and its static extended part. The 
larger FWHM of the precession peaks at I > is likely 
caused by thermally activated transitions between differ- 
ent configurations of F2 stabilized by ST, also resulting 
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FIG. 2: (a) Total microwave power emitted at 1 — 13.5 GHz, 
normalized by the largest possible emitted power, calculated 
as described in the text, (b) FWHM of the first precession 
harmonic vs I. Squares: Py2 at 600 Oe, crosses: Py3 at 
600 Oe, triangles: PyA at 400 Oe. 

in the 1 // noise Q . 

Sample PyA shows precession peaks only at I > 0, 
but no detectable dynamics at / < (Fig. QJe)). The 
dV/dl curve (Fig. QJf)) has a step associated with the 
onset of current-induced dynamics only at I > 0, but 
not at I < 0, suggesting that dynamics is completely 
suppressed at I < 0. It may be possible that a thicker Fi 
has a precession onset outside of the measurement range 
of —8 to 8 mA. However, assuming that 1^ is determined 
predominantly by the dimensions of Fi, one may expect 
that the ratio Iq /^c increases only by about a factor of 
4/3 for PyA compared to Py3, giving Iq w 4 mA well 
within our measurement range. 

Samples with Fi much thicker than F2 exhibited the 
smallest FWHM of precession peaks at H just above H c , 
quickly increasing as H was increased 0. Similarly, the 
400 Oe data for PyA in Fig. [Hg) exhibit smaller 1// noise 
and narrower precession peaks than the 600 Oe data, 
but the smallest FWHM of 165 MHz is still about five 
times larger than in nanopillars with thick Fi. A peak 
in dV/dl data (Fig.QJh)) at / = 6 mA is correlated with 
abrupt termination of precession, which is caused by the 
transition of F2 into a static AP orientation with respect 

to Fr [m. 

The magnetic dynamics can be further characterized 
by comparing the total microwave power emitted by the 
device to the largest possible power emitted by a hy- 
pothetical oscillation between AP and P configurations 
of the magnetic layers, P max = 2I 2 AR 2 R L /(R + R L ) 2 . 
Here, R is the average resistance during oscillation, and 
Rl = 50 f2 is the input impedance of the amplifier. This 
formula was obtained by considering an equivalent elec- 
tric circuit consisting of a dc current source, oscillating 
resistance representing the sample, and a load resistor 
Rl. Fig. HJa) shows the emitted power normalized by 
P max = 10.8/ 2 pW, 27.7I 2 pW, and 29. 7I 2 pW for Py2, 
Py3, and PyA, respectively, where / is in mA. 

For Py2, the maximum relative power P r = P/P ma x — 
0.45 at / = —3 mA corresponds to the in-plane precession 
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FIG. 3: (Color online) (a) PSD for Py2 at H = -600 Oe, 
rotated by 55° with respect to the nanopillar easy axis, (b) 
Integrated power Pi in the first precession harmonic vs its 
peak frequency fi at I < (open symbols) and I > (solid 
symbols). Horizontal bars show FWHM. Arrow on top shows 
the direction of increasing 



angle exceeding 90°, which is remarkable for an extended 
Fi. A significant reduction of P r at / < —5.5 mA is cor- 
related with a dramatic increase of the FWHM of the pre- 
cession peaks (Fig. EJb)). A similar correlation between 
the power and the FWHM of peaks is also observed in 
asymmetric samples with thick Fi, and is likely caused 
by the increasingly complex magnetic dynamics at large 
/ 0. In contrast, sample Pi/3 exhibits simultaneously a 
small P r < 0.06 at all I < (Fig. |5fa)) and a FWHM of 
less than 30 MHz at -5 mA> I > -6.7 mA (Fig. Mb)). 
These behaviors of Py3 are transitional between a large 
power emitted by Py2 and a negligible power emitted 
by PyA at I < 0, suggesting that the dynamics of Fi in 
Py3 is partly suppressed by the same mechanism that 
completely eliminates it in Py4. 

The dependence of the emitted power on tp y at / > is 
opposite to that at / < 0. Sample Py4 shows a maximum 
P r = 0.41 at H = 400 Oe, while Py2 and Py3 show a 
significantly smaller P r < 0.14. The maximum for Py4 
at 600 Oe is P r = 0.28 (not shown). The FWHM for Py4 
first decreases with increasing / < 4.2 mA. Its dramatic 
increase at larger I is correlated with a decrease of P r . 
The FWHM data for Py3 are surprisingly consistent with 
these behaviors, suggesting that the peak broadening is 
governed by the fundamental dynamical properties of the 
bilayers, rather than imperfections of individual devices. 

The dynamics of sample Py2 at / > became more 
coherent when the angle between H and the nanopillar 
easy axis was increased to 55°, as shown in Fig. [3fa). 
This resulted in a decreased amplitude of peaks at / < 0, 
as can be seen from the reduced color scale. We observe a 
similar reduction of the microwave power in all our sam- 
ples exhibiting bipolar precession. The peaks in Fig.[3Ja) 
abruptly broaden at / < —6 mA, and the precession sig- 
nal quickly decreases below the measurement noise floor. 



Figure[3jb) shows that the frequency of precession at this 
point is correlated with the precession onset frequency of 
3.6 GHz at I > 0. For Py3, broadening occurs at larger 
I < — 7 mA (Fig. Elb)), which can be attributed to the 
higher precession frequency of Fi in Py3 than in Pj/2 at 
the same I. We also see a similar broadening of preces- 
sion at I > 5 mA for Py3 and Py4 (Fig. Gfb)), but its 
relation to the dynamical properties of Fi is less clear. 

To understand the implications of our results, we first 
assume that the dynamics of Fi at / > is negligible. ST 
acting on F 2 should result in its precession at I > Iq > 0, 
which is determined only by the polarizing properties of 
Fi and the magnetic moment of F 2 . One can expect a 
modest increase of Iq with decreasing tp y due to the 
weaker spin-polarizing properties of a thin F^ Similarly, 
one may expect that the magnitude of the precession 
onset current Iq for Fi increases with increasing tp y due 
to the larger magnetic moment of Fi . Our results show 
that the precession of each layer rapidly disappears when 
its relative thickness is increased, instead of a gradual 
variation of the onset currents. This dramatic variation 
of the current-induced behaviors cannot be explained in 
terms of independent effects of ST on the two layers, 
and instead requires a self-consistent treatment of the 
simultaneous effects of ST on both layers. 

Analysis of dynamics for magnetic layers with similar 
dynamical properties showed that ST-induced coupling 
can result in precession decoherence explaining the 
broadening of the precession peak at / < when its 
frequency became close to the onset frequency at I > 
(Fig. [3|). Coherent bipolar precession is not possible 
in structures with similar dynamical properties of the 
magnetic layers, because ST-induced coupling completely 
suppresses precession of the thicker layer. 

To address the possibility of bipolar excitations in 
symmetric structures with different dynamical proper- 
ties of layers, we develop a simple analytical model de- 
scribing dynamics of two macrospins Si and S 2 cou- 
pled by ST. Here, the part of the extended layer Fi 
in the point contact area is approximately treated as a 
macrospin Si experiencing an additional effective field 
due to the exchange coupling to the extended film. Lin- 
earizing the Landau-Lifshits-Gilbert-Slonczewski equa- 
tion [12j for small-amplitude precession of unit vectors 
si,2 = Si,2/ |Si,2 1 around the equilibrium orientation 
along the z-axis, we obtain Q 

^77^ = ^1,2^1.2 - r x 2 £l,2 - Ikx 2 (£l - £ 2 ) (1) 

at 

where £i j2 = s\ x _2x + is\y,2y are complex coordinates 
representing the projections of the macrospins on the 
xy plane. The first term on the right describes preces- 
sion around the magnetic field, the second term describes 
relaxation, and the third term describes interaction be- 
tween Si and s 2 due to ST. The coefficients c<Ji. 2 depend 
on £1,2 due to the demagnetizing fields of nanomagnets. 
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By rescaling s x , s y and t, this dependence can be trans- 
ferred to I\2 and kx t 2- Assuming that the precession 
term is dominant, these coefficients can be replaced with 
their averages over the precession cycle. Within these 
approximations, Eqs. ([T]) with constant coefficients can 
be used to analyze the possibility of stable bipolar exci- 
tations in a symmetric bilayer. 

Let's assume that £^2 precess at a common frequency 
w, with a precession phase offset </>, £1 = Aie luJt , and £2 = 
A 2 e lult+i ^ '. Equations ((lj then reduce to two complex 
algebraic relations 

iu = iwi, 2 - r 1>2 T Jfci,a[l - (re # ) ±1 ] (2) 

for the unknowns u>, (f>, r = A2/A1, and I. The value 
of / determined from these relations corresponds to the 
onset of the current-induced dynamics. Assuming for 
simplicity Ti = and subtracting the two equations, we 
obtain cos0 = (fci + k2)/{k\r + k^lr) from the real part, 
and sin0 = {u\ — UJ2) / (I k 2 / r — I kir) from the imaginary 
part. For / > 0, the dynamics of F 2 may be expected 
to dominate, so that r>l. In this case, cos<p « (fci + 
k^j/k^r <C f, and therefore sin0 rs — f for lo\ > u>2, 
giving r = (u>i — u)2)/Iki. From Eqs. lu rs lo 2 , 

and J w T/k 2 - Stability analysis can be performed by 
assuming small deviations Ar and A<f> of the respective 
variables from the stationary values determined above, 
and inserting them into Eqs. {!]). The spectrum of the 
resulting linear differential system lies entirely in the left 
half of the complex plane, confirming the stability of this 
dynamical state. 

From the expressions for r and I above, one can see 
that the condition r ^> 1 is satisfied for similar k\ and k 2 
only if \u)\ — W2I 3> T, i.e. when the dynamical properties 
of the magnetic layers are sufficiently different. For \ui% — 
w 2| < r, the dynamics must instead involve precession of 
both layers with a similar amplitude 0] • Analysis similar 
to that given above also yields a stable dynamical state 
with r <C 1 at / < 0, as long as u\ and L02 are sufficiently 
different. In contrast, assuming r <C 1 at I > or r 3> f 
at / < 0, we obtain a dynamical state which is not stable 
with respect to small perturbations of r and <f>. 

Our analysis supports the possibility of bipolar exci- 
tations in symmetric bilayers with different dynamical 
properties of the layers. The model does not limit the 
relative thicknesses of the two layers, which seems to be 
inconsistent with the rapid suppression of bipolar dy- 
namics in asymmetric samples. We propose two possi- 
ble mechanisms for this suppression. Firstly, small-angle 
precession frequency / w 8 GHz for F2 at H — 600 Oe 
can be calculated using Kittel formula. The onset fre- 
quency of about 4 GHz in our data is rather associated 
with large-amplitude clamshell precession. Small-angle 
elliptical precession is not seen due to the small ampli- 
tude and the large width of the peaks at the precession 
onset (see Fig.[5Jb)). Because of this significant decrease 
of frequency with amplitude, F2 can behave as a highly 



nonlinear oscillator whose frequency can autotune into 
resonance with Fi, suppressing the precession of the lat- 
ter as long as the amplitude of precession of F2 is not lim- 
ited by its larger thickness. A similar argument can be 
used to explain the suppression of precession of F2 by the 
dynamics of Fi. Secondly, for any precession frequency 
of F2 , there is an inhomogeneous dynamical mode of ex- 
tended Fi at the same frequency, which may be critical 
for suppression of the precession of F2 due to the dy- 
namical ST-induced coupling. It is not clear if a similar 
argument can be made for the suppression of the dynam- 
ics of Fi. A more detailed analysis explicitly including 
the nonlinearity of oscillations and inhomogeneous exci- 
tations may be needed to clarify these behaviors. 

In conclusion, we have demonstrated bipolar current- 
induced precession in symmetric magnetic nanopillars 
with different dynamical properties of the two layers. 
Bipolar precession is rapidly suppressed in asymmetric 
devices. This behavior is interpreted in terms of the 
dynamical coupling between magnetic layers induced by 
spin transfer. Our results suggest that coherent pre- 
cession in nanomagnetic devices can be most efficiently 
achieved with significantly different dynamical properties 
of the magnetic layers. Only a small asymmetry between 
their dimensions is then needed to achieve strongly asym- 
metric behaviors of the bilayer. 
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